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Abstract 

Recent studies have thrown up a big enigma. On the one hand 
they point to satellite galaxies rotating faster than they should with 
the usual theory. On the other hand they also go against explana- 
tion of Newtonian gravity with dark matter. In this brief note it is 
argued that the varying G cosmology which explains all the observed 
General Relativistic effects and the Pioneer anomaly reconciles the 
controversial conclusions of the latest observations. 
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1 Introduction 

It is well known that F. Zwicky introduced the concept of dark matter to 
account for the anomalous rotation curves of the galaxies [HE]- The problem 
was that according to the usual Newtonian Dynamics the velocities of the 
stars at the edges of galaxies should fall with distance as in Keplarian orbits, 
roughly according to 

. - M (i) 
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where M is the mass of the galaxy, r the distance from the centre of the 
galaxy of the outlying star and v the tangential velocity of the star. Obser- 
vations however indicated that the velocity curves flatten out, rather than 
follow the law (JT]). This necessitated the introduction of the concept of dark 
matter which would take care of the discrepancy without modifying Newto- 
nian dynamics. However even after nearly eight decades, dark matter has not 
been detected, even though there have been any number of candidates pro- 
posed for this, for example SUSY particles, massive neutrinos, undetectable 
brown dwarf stars, even black holes and so on. 

Very recent developments are even more startling. These concern the ro- 
tating dwarf galaxies, which are satellites of the Milky Way [3J H]. These 
studies throw up a big puzzle. On the one hand these dwarf satellites can- 
not contain any dark matter and on the other hand the stars in the satellite 
galaxies are observed to be moving much faster than predicted by Newtonian 
dynamics, exactly as in the case of the galaxies themselves. Metz, Kroupa, 
Theis, Hensler and Jerjen conclude that the only explanation lies in rejecting 
dark matter and Newtonian gravitation. Indeed a well known Astrophysicist, 
R. Sanders from the University of Groningen commenting on these studies 
notes [5], "The authors of this paper make a strong argument. Their result 
is entirely consistent with the expectations of modified Newtonian dynam- 
ics (MOND), but completely opposite to the predictions of the dark matter 
hypothesis. Rarely is an observational test so definite." Moreover Vahe Pet- 
rosian of the Kavli Institute in the February 10 issue of the Astrophysical 
Journal, rules out dark matter on the basis of studying the interaction of 
electrons at the galactic edge with starlight. 

In this note we point out that this could indeed be so, though not via Mil- 
grom's ad hoc modified dynamics according to which a test particle 
at a distance r from a large mass M is subject to the acceleration a given by 

a 2 /a = MGr' 2 , (2) 

where ao is an acceleration such that standard Newtonian dynamics is a 
good approximation only for accelerations much larger than ao- The above 
equation however would be true when a is much less than ao- Both the 
statements in (jSJ) can be combined in the heuristic relation 

fi(a/a )a = MGr' 2 (3) 

In (j3J) fi(x) ~ 1 when x » 1, and/i(x) ~ x when x << 1. It must be 
stressed that (J2J) or are not deduced from any theory, but rather are an ad 
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hoc prescription to explain observations. Interestingly it must be mentioned 
that most of the implications of Modified Newtonian Dynamics or MOND 
do not depend strongly on the exact form of \i. 

It can then be shown that the problem of galactic velocities is now solved 
[7J [U [9] [10] . Nevertheless, most physicists are not comfortable with MOND 
because of the ad hoc nature of (|2J) and 

2 Varying G Dynamics 

We now come to the cosmological model described by the author in 1997 
(Cf.ref.[TTJ [2] and several references therein), in which the universe, under 
the influence of dark energy would be accelerating with a small accelera- 
tion. Several other astrophysical relations, some of them hitherto inexpli- 
cable such as the Weinberg formula giving the pion mass in terms of the 
Hubble constant were also deduced in this model (Cf.also ref.[12] and refer- 
ences therein). While all this was exactly opposite to the then established 
theory, it is well known that the picture was observationally confirmed soon 
thereafter through the work of Perlmutter and others (Cf.ref. [12]). Inter- 
estingly, in this model Newton's gravitational constant varied inversely with 
time. 

Cosmologies with time varying G have been considered in the past, for ex- 
ample in the Brans-Dicke theory or in the Dirac large number theory or by 
Hoyle [13l [HI EE [161 E] • m the case of the Dirac cosmology, the motiva- 
tion was Dirac's observation that the supposedly large number coincidences 
involving N ~ 10 80 , the number of elementary particles in the universe had 
an underlying message if it is recognized that 



where T is the age of the universe. Equation (j4j) too leads to a G decreasing 
inversely with time as we will now show. We follow a route slightly different 
from that of Dirac. 
From (H]) it can easily be seen that 



where r is a typical Compton time of an elementary particle ~ 10 sees, 
because T, the present age of the universe is ~ 10 17 secs. We also use the 



y/N oc T 



(4) 



T = VNt 



(5) 
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following relation for a uniformly expanding Friedman Universe 

^2 = ^ GjR 2 p (6) 

where R is the radius of the universe and p its density. We remember that 

p = and M = Nm (7) 

AtcR 3 

where M is the mass of the universe, and m is the mass of an elementary 
particle ~ 10~ 25 gm (Cf.ref.fT5]). 

Use of (H) in (16]) leads to another well known relation [19] 



c 2 



because R = c. Further dividing both sides of (jSJ) by c we get the famous 
Weyl-Eddington relation 

R = y/Nl (9) 

where I — r/c is a typical Compton length ~ 10~ 13 cms. 
Use of (ED and © in © now leads to 

Equation (|T0|) gives the above stated inverse dependence of the gravitational 
constant G on time, which Dirac obtained. On the other hand this same 
relation was obtained by a different route in the author's dark energy - fluc- 
tuations cosmology in 1997. This work, particularly in the context of the 
Planck scale has been there for many years in the literature (Cf.[T2l [2] ED] 
and references therein). Suffice to say that all the supposedly so called ac- 
cidental Large Number Relations like (Q as also the inexplicable Weinberg 
formula which relates the Hubble constant to the mass of a pion, follow as 
deductions in this cosmology. The above references give a comprehensive 
picture. 

The Brans-Dicke cosmology arose from the work of Jordan who was moti- 
vated by Dirac's ideas to try and modify General Relativity suitably. In this 
scheme the variation of G could be obtained from a scalar field <fi which would 
satisfy a conservation law. This scalar tensor gravity theory was further de- 
veloped by Brans and Dicke, in which G was inversely proportional to the 
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variable field <fi. (It may be mentioned that more recently the ideas of Brans 
and Dicke have been further generalized.) 

In the Hoyle-Narlikar steady state model, it was assumed that in the Machian 
sense the inertia of a particle originates from the rest of the matter present 
in the universe. This again leads to a variable G. The above references give 
further details of these various schemes and their shortcomings which have 
lead to their falling out of favour. 

In any case, our starting point is, equation ( TTOj) where T is time (the age of 
the universe) and Go is a constant. Furthermore, other routine effects like 
the precession of the perihelion of Mercury and the bending of light and so 
on are also explained with ( II Op as will be briefly discussed below. We will 
also see that there is observational evidence for ([TO]) (Cf. also [21] which 
described various observational evidences for the variation of G, for exam- 
ple from solar system observations, from cosmological observations and even 
from the palaeontological studies point of view). 

With this background, we now mention some further tests for equation ([TO]) . 
This could explain the other General Relativistic effects like the shortening 
of the period of binary pulsars and so on (Cf.ref. [121 [21 1221 123] and other 
references therein). Moreover, we could now also explain, the otherwise in- 
explicable anomalous acceleration of the Pioneer space crafts (Cf.ref.[2] for 
details). We will briefly revisit some of these effects later. 
We now come to the problem of galactic rotational curves mentioned earlier 
(cf.ref.pQ). We would expect, on the basis of straightforward dynamics that 
the rotational velocities at the edges of galaxies would fall off according to 



which is (JTJ) . However it is found that the velocities tend to a constant value, 



This, as noted, has lead to the postulation of the as yet undetected additional 
matter alluded to, the so called dark matter. (However for an alternative view 
point Cf. [23]). We observe that ( ITU]) can be written for an increase t,in time, 
small compared to the age of the universe, now written as to 



v ~ 300km/ sec 



(12) 




(13) 
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Using ( 1T5|) . let us consider the gravitational potential energy V between two 
masses, vn\ and m 2 by: 

_ Gm 1 m 2 _ Go mim 2 
r to r 



After a time £ this would be, by ([[3 

Go A 

to V to 



y = ^fi-^^ (15) 



Equating (TH|) and ([15]) we get, 



- 1 vFT ) < 16) 



The relation ( |T6l) has been deduced by a different route by Narlikar pQ . From 
(|T6|) it easily follows that, 

a = (f G - f) « i(ir + 2r ) « -2 J (17) 

to t Q 

as we are considering intervals t small compared to the age of the universe and 
nearly circular orbits. In ( TT7|) . a or the left side of (|T7|) gives the new extra 
effect due to (|T3l) and (|T6|) . this being a departure from the usual Newtonian 
gravitation. Equation (fT7|) shows (Cf.ref[22] also) that there is an anomalous 
inward acceleration, as if there is an extra attractive force, or an additional 
central mass. 

So, introducing the extra acceleration f[T7|) . we get, 



r 2 t 

From ( fT8|) it follows that 



GMm 2mr mv 2 , . 

~^ + — ^— ( 18 ) 



'2r 2 GM\ 1/2 



So ( TL9l) replaces ([I]) in this model. This shows that as long as 

2r 2 GM 



(19) 



« , (20) 



tl r 
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Newtonian dynamics holds. But when the first term on the left side of ( T2"U]) 
becomes of the order of the second (or greater), the new dynamical effects 
come in. 

For example from f fl9|) it is easily seen that at distances well within the edge 
of a typical galaxy, that is r < 10 23 cms the usual equation ( ITT]) holds but as 
we reach the edge and beyond, that is for r > 10 24 cms we have v ~ 10 7 cms 
per second, in agreement with f[T2"]) . In fact as can be seen from ([!§]) . the 
first term in the square root has an extra contribution (due to the varying 
G) which exceeds the second term as we approach the galactic edge, as if 
there is an extra mass, that much more. 

We would like to stress that the same conclusions will apply to the latest 
observations of the satellite galaxies (without requiring any dark matter). 
Let us for example consider the Megallanic clouds [25]. In this case, as we 
approach their edges, the first term within the square root on the right side 
of (]T9]) or the left term of (120]) already becomes of the order of the second 
term, leading to the new non Newtonian effects. 



3 Remarks 

We have already noted that the varying G dynamics given in (TTO]) explains 
all the General Relativistic effects. This work has been available in the liter- 
ature for many years (Cf.[22][I2][2ni[2] for full details and further references). 
However we repeat some examples briefly to give an idea. 
We could explain the correct gravitational bending of light. Infact in Newto- 
nian theory too we obtain the bending of light, though the amount is half that 
predicted by General Relativity [JJ |2S1 [23 I2B]- In the Newtonian theory we 
can obtain the bending from the well known orbital equations (Cf.also[2l?]). 

- = ^-{l + ecosQ) (21) 
r L z 

where M is the mass of the central object, L is the angular momentum 
per unit mass, which in our case is 6c, h being the impact parameter or 
minimum approach distance of light to the object, and e the eccentricity of 
the trajectory is given by 

C = 1 + G^ < 22 > 
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For the deflection of light a, if we substitute r = ±00, and then use (122]) we 

get 

2GM , , 

<* = 13" (23) 

This is half the General Relativistic value. 

We now note that the effect of time variation of r is given by equation 
( I16|) (cf.ref. |22j). Using this the well known equation for the trajectory is 
given by, 

GM t „/t\ 2 

where u = \ and primes denote differenciation with respect to 0. 
The first term on the right hand side represents the Newtonian contribution 
while the remaining terms are the contributions due to ( ITOl) . The solution of 
( |24|) is given by 



GM 

u = ~u~ 



1 -'r (COS {{1 —](-) + 

2t 



(25) 



where u is a constant of integration. Corresponding to — 00 < r < 00 in the 
Newtonian case we have in the present case, —to < t < to, where to is large 
and infinite for practical purposes. Accordingly the analogue of the reception 
of light for the observer, viz., r = +00 in the Newtonian case is obtained by 
taking t = to in ( 1251) which gives 

GM /6 \ 

u = —j-p- + ecos I — + oj I (26) 

Let us compare ( 1261) with the Newtonian solution. This shows that the 
Newtonian 6 is replaced by ®, whence the deflection obtained by equating 
the left side of (|26|) to zero, is 



where e is given by f[22|) . The value of the deflection from fl27j) is twice the 
Newtonian deflection given by fl23|) . That is the deflection a is now given 
not by ( 1231) but by the formula, 

AGM , . 

= - w-- < 28 ' 
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The relation (128 j) is the correct observed value and is the same as the General 
Relativistic formula which however is obtained by a different route [2~g| I50| 



Next, we come to the inexplicable anomalous accelerations of the Pioneer 
spacecrafts already alluded to, which have been observed by J.D. Anderson 
and coworkers at the Jet Propulsion Laboratory for well over a decade [521155] 
and have posed a puzzle. This can be explained in a simple way as follows: 
In fact from the usual orbital equations we have 

CM ... CM 

2V' 



vv ~ —7^—7(1 + ecosQ) 2~^(1 ecos ®) 



v being the velocity of the spacecraft and t is the time in general. It must be 
observed that the first term on the right side is the new effect due to ({TO]) . 
There is now an anomalous acceleration given by 

a r = Wanom = -77; 1 1 + ecosB) 

2trv 

1 + e) 3 



2t\ 
where 

A = r 4 6 2 

If we insert the values for the Pioneer spacecrafts we get 

a r ~ — 10~ 7 cm/sec 2 

This is the anomalous acceleration reported by Anderson and co-workers. 
We will next deduce that equations like (|T6|) explains correctly the observed 
decrease in the orbital period of the binary pulsar PSR 1913 + 16, which has 
also been attributed to as yet undetected gravitational waves [55] . 
It may be observed that the energy E of two masses M and m in gravitational 
interaction at a distance L is given by 

rn GMm , . 

E = = constant (29) 

L 

We note that if this energy decreases by any mechanism, for example by 
the emission of gravitational waves, or by the decrease of G, then because 
of (|29p . there is a compensation by the decrease in the orbital length and 
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orbital period. This is the standard General Relativistic explanation for the 
binary pulsar PSR 1913 + 16. We will show that the same holds good, if we 
are given instead, ffTUj) . That is, we will not invoke gravitational waves. In 
this case we have, from f l29|) 

— = — : — — = const. (30) 

LJ LJ 

We can now write, for a time increase t, 

ji = Hq — Kt (31) 

where we have 

K = fi (32) 

In ( 132]) fi can be taken to be a constant in view of the fact that G varies very 
slowly with time. Specifically we have 

G(T + t)= G(T) -*£ + ££ + ...« G(T) - t% (33) 

where T is the age of the universe and t is an incremental time. Whence 
using (133]) . K in fl32|) is given by 

Ko,- 
T 

and so 

C 

rp2 



So ( 1301) requires 

L = L (l - aK) 

Whence on using ([3*1]) we get 

a = — (34) 

Let us now consider t to be the period of revolution in the case of the binary 
pulsar. Using (151]) it follows that 

SL = (35) 
Mo 
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We also know (Cf.ref. [33]) 



271 , 9 27T , . 

t = — L = 36 

h y/Ji ^ ^ 

i 2 = , 37 

/i being the usual unit angular momentum and /i has the units gmcm^sec^ 1 . 
Using ( 1351) . ( 1361) and ( 1371) . a little manipulation gives 

2t 2 K 

5t=-— (38) 

(135j) and (1581 show that there is a decrease in the size of the orbit, as also in 
the orbital period. Such a decrease in the orbital period has been observed 
in the case of binary pulsars in general (351 EE] ■ 

Let us now apply the above considerations to the specific case of the binary 
pulsar PS R 1913 + 16 observed by Taylor and coworkers (Cf.ref. (36]). In this 
case it is known that, t is 8 hours while v, the orbital speed is 3 x 10 7 cms 
per second. It is easy to calculate from the above 

Ho = 10 4 x v 3 ~ 10 26 

which gives M ~ 10 33 gms, which of course agrees with observation. Further 
we get using (j4j) and ( |3T!) 

At = r]X 10~ 5 sec/yr, r] <^ 8 (39) 

Indeed (139]) is in good agreement with the carefully observed value of r\ ~ 7. 5 
(Cf.refs. [251136]). 

It should also be remarked that in the case of gravitational radiation, there 

are some objections relevant to the calculation (Cf.ref. [3 5J). 

Finally, we may point out that a similar shrinking in size with time can 

be expected of galaxies themselves, and in general, gravitationally bound 

systems. 

To consider the above result in a more general context, we come back to the 
well known orbital equation (3^ 



d 2 u/dQ 2 + u = fi /h 2 (40) 
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where /xo = GM and u is the usual inverse of radial distance. 

M is the mass of the central object and h = r 2 d<d/dt - a constant. The 

solution of ( I40p is well known, 

lu = 1 + ecosQ 

where I = h 2 / hq. 

It must be mentioned that in the above purely classical analysis, there is no 
precession of the perihelion. 

We now replace /x by /x and also assume /x to be varying slowly because 
G itself varies slowly and uniformly, as noted earlier: 

(i — dfi/dt — K,a constant (41) 

remembering that K ~ 0(1/T 2 ) and so can be neglected. 

Using (T4TT) in fT40l) and solving the orbital equation (T40I) . the solution can 
now be obtained as 

u=l/l+ (e/l)cos& + Kl 2 Q/h 3 + Kl 2 eQcos&/h 3 (42) 

Keeping terms up to the power of V and (K/fi ) 2 , the time period V for 
one revolution is given to this order of approximation by 

r = 2nL 2 /h (43) 

From (JH) 

L = l-i 

h 3 



L = l- —— (44) 



Substituting (jSJ) in ( T431 we have 

The second term in fT4"5l) represents the change in time period for one revo- 
lution. The decrease of time period is given by 

St = 8tv 2 1 3 K/ijI (46) 

The second term in (jUJ) indicates the decrease in latus-rectum. 
For one revolution the change of latus-rectum is given by 

51 = 2nKl 2 - 5 / f i 1 5 (47) 
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In the solar system, we have, 



K = 898800 cm gm 

Using K and /xq to find the change in time period and the latus rectum in 
the varying G case by substituting in (146]) and (|47|) respectively for Mercury 
we get 

ST = 1.37 x 10- 5 sec/rev 

SI = AMcm/rev (48) 

We observe that the equations (146}) . ( 14T]) or (148]) show a decrease in distance 
and in the time of revolution. If we use for the planetary motion, the General 
Relativistic analogue of (HI)]) , viz., 

— +^-(1 + 3/^), 

then while we recover the precession of the perihelion of Mercury, for example, 
there is no effect similar to (146]) , ( 147]) or (148]) . On the other hand this effect is 
very minute- just a few centimeters per year in the case of the earth- and only 
protracted careful observations can detect it. Moreover these changes could 
also be masked at least partly, by gravitational and other perturbations. 
However as noted, the decrease of the period in (T46]) has been observed in 
the case of Binary Pulsars. 



4 Conclusion 

The point is that the above varying G scheme described in ( flOl) or ( [19]) re- 
produces all the effects of General Relativity as noted above, as also the 
anomalous acceleration of the Pioneer space crafts in addition to the conclu- 
sions of MOND regarding an alternative for dark matter, and is applicable 
to the latest observations of satellite galaxies. The satellite galaxy rotation 
puzzle is thus resolved. 
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